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0\ ■ 1. INTRODUCTION 
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"Y* 'i Albert Einstein found some radiative solutions to the field equations of gen- 

^ ' eral relativity and deduced the existence of gravitational waves, which could be 

interpreted as ripples in the curvature of spacetime propagating with the speed 
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of light. In a subsequent paper he computed the total power emitted in the 
form of gravitational waves by an isolated "Newtonian" source, and found that 
; ^ ' this power depends quadratically on the variations of the quadrupole moment 

■ of the source in a celebrated formula known today as Einstein's quadrupole for- 

mula. The fact that gravitational waves from an isolated source are dominantly 
quadrupolar is a consequence of the equivalence principle, which implies the 
conservation (or linear variation in time) of the source's monopole and dipole 
moments as the consequence of the equations of motion. 

Existence of gravitational radiation was confirmed by Eddington who 
showed that the propagation at the speed of light of weak gravitational waves 
on a flat background has an intrinsic meaning (independent of the choice of a 
coordinate system), and by Bondi Q who proved that it is possible to construct 
a detector of gravitational waves which is heated when the wave passes through 
it. Surely gravitational radiation thus carries energy. Furthermore this energy 
is extracted from the mass-energy of the source (as measured from infinity), 
which has been shown by Bonnor ^ and Bondi, van der Burg and Metzner Q 
(see also Sachs [0) to be a monotonically decreasing function of time during the 
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emission. Before emission this mass-energy reduces to the mass of Arnowitt, 
Deser and Misner Q which represents the total mass-energy of the source plus 
the contribution of the waves to be emitted. Some theorems due notably to 
Shoen and Yau |l^ have established the positivity of both the ADM and 
Bondi mass of an isolated system. This proves in particular that the mass of 
the source bounds from above the total amount of energy which it can emit in 
the form of radiation. All these theoretical works show that the classic vision 
of a wave carrying energy, which is extracted from its source and can be put 
down an antenna, is legitimate in the case of gravitational radiation. 

On the observational side gravitational radiation is known to exist thanks to 
the discovery by Hulse and Taylor of the binary pulsar PSR 1913-1-16. Five 
years of accurate timing observations of this pulsar have yielded the conclusion 
(by Taylor, Fowler and Mc CuUoch |l^) that the orbital period of the binary 
system of the pulsar and its companion is steadily decreasing, and therefore that 
by Kepler's third law the system is loosing energy. There is exact agreement 
with the value of the flux of energy by gravitational radiation as predicted by 
the Einstein quadrupole formula ||l^-|T^. This effect can also be seen as due 
to the reaction forces acting on the orbit in reaction to the emission of gravita- 
tional waves. Radiation reaction forces in general relativity and the associated 
energy balance equation have been computed by Chandrasekhar and Esposito 
(see also |l^-||l^). Burke and Thorne and many others authors 

pG|. None of these computations, however, gave a complete description of 
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the dynamics of gravitating systems, nor were applicable to systems contain- 
ing strongly self-gravitating bodies such as neutron stars. The first complete 
derivation of the dynamics of the binary pulsar, up to the level where radiation 
reaction effects appear, was obtained by Damour and Deruelle l^-fs^ (see 
the lecture of T. Damour in this volume). 

Experimental research on gravitational radiation is presently very active with 
the development of new technologies for bar detectors, and most importantly 
with the construction of two large scale laser interferometric detectors for ob- 
servation of the waves in the frequency bandwidth between w 10 Hz and w 1000 
Hz : the American LIGO detector [B6l, and the Franco-Italian VIRGO detector 



37 . These detectors should observe at least one type of source, the inspiralling 
compact binary. This source is composed of two neutron stars or black holes 
spiralling very rapidly around each other in the last rotations preceding their 
final coalescence. Orbital velocities are much larger than in the binary pulsar 
for instance, which will coalesce with its companion in few hundreds of millions 
years. 

For very relativistic sources like inspiralling compact binaries, the precision 
given by the Einstein quadrupole formula for the energy in the waves is un- 
sufficient. Similarly unsufhcient is the precision of the radiation reaction force 
of Burke and Thorne for instance. What is required is a relativistic or post- 
Newtonian formalism (involving an expansion when the speed of light goes to 
infinity) for both the emission and reaction of waves from isolated sources with 
substantially large internal velocities. However as we do not know what are the 
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sources which will actually be detected by VIRGO and LIGO, it is important 
to have a formalism which is sufficiently general (not limited, for instance, to 
inspiralling compact binaries). 

Pionneering work for extending the Einstein quadrupole formula (and for- 
mulas alike) to post-Newtonian order is due to Epstein and Wagoner [3^ . Note 
however that the post-Newtonian expansion is limited to the near-zone of the 
source (see Section 2.1). Thus, for studying the radiation field one needs a pri- 
ori to supplement the post-Newtonian expansion by some another method. The 
post-Minkowskian method (expansion when the gravitational constant G goes 
to zero) is adequate for that purpose. Pionneering work on post-Minkowskian 
expansions is due to Bertotti and Plebanski [§9 L Bonnor and collaborators 
and Thorne and collaborators |44f@. Notably Bonnor § in- 
troduced the idea of considering the post-Minkowskian expansion outside the 
source where a multipolar decomposition of the field simplifies the resolution 
of the field equations. 

Our aim in this lecture is to discuss, with emphasis on recent developments, 
a particular post-Newtonian (and post-Minkowskian) formalism for the gen- 
eration of gravitational waves by a relativistic source, and the reaction of the 
waves on the source. Part of this formalism is issued from the line of research 
initiated by Bonnor and collaborators ||], |40 Initial work was inspired 

by Thome's review on multipole expansions ||46|. The formalism is sufficiently 
mature to tackle the problem of inspiralling compact binaries, and we present 
at the end of the lecture the results which have been obtained so far in this ap- 
plication. The work on this post-Newtonian formalism started with the present 
author's PhD thesis (advised by T. Damour and published in |4^). The list of 
references is provided at the end of the lecture |4^~[0. A very useful auxiliary 
technical development was brought by T. Damour and B. Iyer |^8| . Applica- 
tions to binary systems (and, most importantly, inspiralling compact binaries) 
were done in collaboration with G. Schafer ]59|, 60t|, T. Damour and B. Iyer 
]6l]] , and partly with C. Will and A. Wiseman [|62, |63t| who used also their own 
method based on the Epstein and Wagoner approach jsl, |63]. A. Wiseman 



pq, 66 1 and L. Kidder et al pTl B8| made also applications of the formalism. 



The present lecture is strongly biased by the author's point of view, so we 
would like to refer here to other articles for different approaches and methods, 
or simply for different points of views. See in particular the Les Houches book 
of 1982 on Gravitational Radiation (chapters by Thorne, Damour, Walker, 
Choquet-Bruhat, Friedrich, Ashtekar), the chapters of Thorne and Damour 
in the book on 300 years of gravitation (see in particular the discussions 
on approximation methods), the review by Damour in the Cargese book on 
Gravitation in Astrophysics (see especially the discussion on the various 
quadrupole laws and equations), the chapters by Iyer and Schmidt in Advances 
in Gravitation and Cosmology |7^, several reviews by Thorne |7^-[^, Schutz 
1?^ , [ttI and Bonazzola and Marck |Q more devoted to detectors and the 
sources of waves, and one review by Will ||7^ on the post-Newtonian approach 
to gravitational radiation. See also the lectures of J. Bicak and T. Damour in 



4 



L. Blanchet 



this volume. 

Wc will proceed didactically starting in Section II with several useful defi- 
nitions of regions around an isolated source, some recalls from the quadrupole 
formalism, and a general outline of the method. Sections III and IV are devoted 
respectively to some relevant properties of radiative gravitational fields and the 
generation of waves by the source. Section V deals with the radiation reaction 
effects occuring within the source. Finally Section VI applies the results of this 
wave generation formalism to inspiralling compact binaries. 

2. BASICS AND OUTLINE OF THE METHOD 
2.1. Regions eiround an isolated source 

All over this lecture, we consider a general isolated source of gravity described 
by a matter stress-energy tensor with components T^"{x, t) in some coordinate 
system (x, i) covering the source (Greek indices take the values 0,1,2,3, and 
Latin 1, 2, 3). Let a be the radius of a sphere which totally encloses the source. 

We shall make the restriction that this isolated source is slowly moving, 
in the sense that the typical internal velocity inside the source divided by 
the speed of light c, i.e. the typical ratio 7"°*/!"'^", defines a small parameter 
e pa 1/c referred to as the slow motion or post-Newtonian parameter. We 
assume further that the typical stresses inside the source divided by the energy 
density in the source, i.e. the typical ratio T^J /T^o^ jg q£ order of the square of 
the parameter s. Finally we consider a purely gravitational problem, so that 
the source is self-gravitating (its internal motions are driven by gravitational 
forces). In this case the ratio GM/c^a, where M is the total mass-energy (or 
ADM mass) of the source, is also of order e^. This implies that the gravitational 
field is weak (and of order e^) everywhere inside and outside the source. 

Note that the parameter e is required to be small, but not exceedingly small. 
For instance it could be as large as say 0.2 or 0.3. This is in order to allow for 
the very interesting case of inspiralling compact binaries where s can reach 0.3 
during the last rotations. Thus we are here considering truly relativistic sources, 
whose accurate description necessitates controlling many relativistic corrections 
(how many depending on the exact range of magnitude of e). However, what 
we do not consider in this lecture is the case of a fully relativistic source, which 
is made of massless particles moving at the speed of light, or, to give another 
example, a cosmic string whose internal stress is of the same order as its energy 
density. 

A useful characterization of slowly moving sources is that their spatial ex- 
tension is small (and of order e) as compared to one typical wavelength of the 
radiation. Indeed the wavelength is given by A = cP where P is a period of 
motion in the source. But a ~ vP with v ~ sc, thus cP a/s and therefore 
a/X £ which is indeed the statement above. Note that since GM/a(? k, 
we have also GM/\c^ k. e^. 
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It is convenient for slowly moving sources to introduce an interior domain 
Di, sometimes called also the near zone, defined by Di = {(x, t), |x| < r^}, 
where the radius is such that ~ eX and ri > a. The near zone Di is small 
with respect to the wavelength of the radiation and covers entirely the source. 
This is possible only for a slowly moving source. We choose to be strictly 
larger than a (instead of being a itself) for later convenience. 

The near-zone is the domain where one can confidently use the post-Newton- 
ian expansion. The real precision of some post-Newtonian expression will be 
exactly given, in Di, by the formal order in e of the neglected terms. Note 
that all powers of 1/c must be taken into account in finding the magnitude of 
a term in terms of e ~ 1/c, including the ones which arise from the temporal 
gradients do = c~^d/dt (which are really of order e with respect to the spatial 
gradients di). This is clear from the definition of the near zone, where the field 
is quasi-static and propagation effects are small. These considerations are also 
familiar from electromagnetism. 

Having clarified the concept of near-zone Di, it is now necessary to introduce 
an exterior domain defined by Df, = {(x, i), |x| > r^}, where the radius is 
chosen to be strictly between a and r^, i.e. a < Vf, < ri. This choice, which 
can always be done for slowly-moving sources, is in order that the intersection 
between Di and D^ (the exterior part of the near zone Di n D^,) exists. 

Included in the exterior domain D^, we also consider the so-called exponen- 
tially far wave zone, namely the domain D^ = {(x, t), |x| > r^^}, where r^^ is a 
radius such that « Ae'^'^ /GM ^ jy^^ where the observer will be located, 
one can expand the field in powers of the distance of the source, and keep the 
leading order term in this expansion. The error done in assuming this will be 
negligible for future astrophysical sources of radiation. 

As is well-known from the study of the Schwarzschild solution, the actual 
space-time light cones at large distances from the black hole deviate by a log- 
arithmic term from the "flat" light cones associated with some coordinate sys- 
tems more suitable near the source (for instance a harmonic coordinate system, 
or a Schwarzschild- like one) . The definition of D^, as being the exponentially 
far wave zone is chosen precisely in order that this effect becomes apprecia- 
ble there. Our exponentially far zone corresponds to the distant wave zone in 
the terminology of Thorne Q| , who introduces also a local wave zone starting 
several wavelength apart from the source. We shall not need here to make 
this distinction. It will be sufficient to consider a minimal splitting into three 
regions Di, D^ and D^. See Thorne [Q for physical discussions of the effects 
which arise in different regions, and for a more general situation of a possibly 
strong-field source and a non-flat gravitational background. 

Let us comment that all the radius a, ri, and which have been in- 
troduced will in fact never appear explicitly in the results presented below. 
We shall only use the existence of the various domains as determined by these 
radius, for instance the fact that the intersection Di n De is not empty. This 
will permit us to prove some statements about the construction of solutions, 
for instance that one can licitly match the outer field to the inner field of the 
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source in DinDe, but the precise locations of the radius are unimportant (these 
radii are anyway loosely defined). 

2.2. The quadrupole formalism 

By quadrupole formalism we mean the lowest order formalism, in the post- 
Newtonian expansion e ^ 0, for the generation of gravitational radiation from 
the source, and also for the reaction of the radiation onto the source. Although 
this may seem to be a little paradoxical, the quadrupole formalism can thus be 
viewed as a Newtonian formalism. Relativistic corrections to the formalism are 
referred to as post-Newtonian corrections accordingly to which post- Newtonian 
order in the equations of motion of the source (beyond the usual Newtonian 
acceleration) is needed in order to reduce all accelerations with consistent accu- 
racy. For instance, the first-post-Newtonian (in short IPN) formalism retains 
all terms in the radiation field and in the reaction which can be computed 
consistently using the IPN equations of motion, which include the corrections 

beyond the Newtonian force (using the standard practice that the nth post- 
Newtonian order refers to the order e^"). 

The quadrupole formalism for the wave generation expresses the gravita- 
tional field hjj^ = {gij — 6ij)'^'^ (where g.ij denotes the spatial part of the 
covariant metric and Sij the Kronecker symbol) in a transverse and traceless 
(TT) coordinate system as 

^r-^^^.-(N){^(J^-i?/c) + 0(s)}+0(i,) . (1) 

The coordinate system (X,T) is centered on the source, with R = |X| the 
distance to the source and N — X/i? the unit direction from the source to the 
observer (later we shall need to be more precise about the exact definition of 
this coordinate system). The retarded time at the observer position in is 
T — R/c, and terms of order in the distance of the source are neglected 

(in addition to the post-Newtonian terms of order e). In front of (1) appears 

V^jkrnCN) = iS,k - N,Nk){S,m ~ N^N,-,) - i(5,, - N,N,)iSkrn - NkN^) (2) 

which is the TT projection operator onto the plane orthogonal to the direction 
N. The quadrupole moment of the source takes the familiar Newtonian form 

Q,,{t) = J d^x p(x,<) i^x'x' - y^'^^^ , (3) 

where p denotes the Newtonian mass density of the source (in some coordinate 
system which we shall also need to specify later). The quadrupole moment (3) 
is taken to be tracefree {SijQij = 0). This is not important in the waveform (1) 
because at the quadrupolar level the TT operator in front cancels any possible 
trace in the moment, but it will be very convenient in higher orders to choose 
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all multipole moments to be symmetric and tracefree (STF) with respect to all 
their indices. 

By differentiating with respect to time the waveform (|l|)-(|^), squaring the 
result using standard formulas and integrating over all the directions N, one 
obtains the total power in the gravitational waves emitted by the source, 

r - — / <^^Qij d^Qtj ri(r'^\ 

This result is the famous Einstein quadrupole formula Q , which neglects post- 
Newtonian terms of order e^. It was derived originally by Einstein under the 
very restrictive assumption that the motion of the source is non gravitational. 
Then Landau and Lifschitz ^ showed that the quadrupole formula ap- 
plies also to a gravitationally bound system, for instance a Newtonian binary 
star system. Many improvements in the derivation of this formula have been 
obtained later, notably by Fock who made use of the idea of matching 
between the inner field of the source and some exterior field in Di D D^- See 
for a completely satisfactory derivation of the quadrupole formula (and its 
first relativistic corrections) using Fock's idea. The notation C in (Q) stands 
for the total luminosity of the source in gravitational waves, by analogy with 
the total luminosity of a star in electromagnetic waves. 

It is straightforward using the expression of the quadrupole moment of the 
source to rewrite the quadrupole formula (Q) in the form of a balance equation 
relating the mechanical energy loss in the source to the total flux of energy in 
the waves. Indeed we can easily rewrite (jj) in the form 

The superscript (n) denotes n time derivatives. The second term in the left- 
hand-side of (H), which has been added to the energy flux, is a total time 
derivative and therefore its time average is numerically small in the case of 
quasi-periodic motion (see e.g. for a discussion). On the other hand, we 
can interpret the right-hand-side as the power extracted within the source by 
the radiation reaction force density ^^^nc^ where T^.^^^ is given explicitly by 

•^;cac = P 9.Kcac+0(£') ; ^reac(x, = " ^ {^'' Q^f (t) + 0{e^)} . (6) 

This expression for the radiation reaction force represents the gravitational 
analogue of the damping force of electromagnetism, and was obtained by Burke 
and Thorne |2^-|Q. Note that gravitational radiation reaction is inherently 
gauge-dependent, so the expression of the force depends on the coordinate 
system which is used (here, a Burke-Thorne coordinate system whose definition 
can be found for instance in Chapter 36 of Misner, Thorne and Wheeler [^). 
Chandrasekhar and Esposito have independently obtained another (and 
more intricated) expression of the reaction force, and Miller [Q has shown 
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how one can transform one expression into the other by means of a coordinate 
transformation (see Schafer [ p5| for a comparison between various expressions 
of the reaction force). 

When deahng with relativistic sources having e ^ 0.3 (say), the post-Newton- 
ian terms which are neglected in the quadrupole formulas both for the wave 
generation (0)-(^ and radiation reaction (||)-(||) will become important. Our 
aim in this lecture is to compute these terms. 

2.3. Outline of the method 

Generalizing the previous quadrupole formalism for application to relativistic 
sources entails solving the mathematical problem of the resolution of the Ein- 
stein field equations with due taking into account of all the nonlinear effects 
entering some high-order post-Newtonian approximation we are considering. 

The general method we use can be decomposed into several different steps. 
As a first step one investigates the gravitational field generated by the isolated 
source in all the exterior domain (which comprises the far wave zone D^). 
The field in Dg is a solution of the vacuum field equations and is computed 
using the combination of two methods : 

(i) a formal infinite post-Minkowskian approximation method for the metric 
field (or perturbation around Minkowski's space-time, or also expansion in 
powers of Newton's constant G) — this method is valid wherever the field is 
weak ; 

(ii) a decomposition of each coefficients of the latter post-Minkowskian ex- 
pansion into multipolc moments (or simply expansion in spherical harmonics), 
which is valid only in the exterior Ug. 

The precise assumptions underlying the methods (i)-(ii) can be found in Sec- 
tion 1.2 of j4^. An iterative resolution {d la Bonnor jsj) of the field equations 
permits us in principle to control, and if necessary to compute, all the non- 
linearities of the exterior field. Spherical harmonics are expressed using the 
symmetric and tracefree (STF) notation which is especially convenient when 
dealing with the decomposition of tensors. This leads also to transparent defi- 
nitions of the multipole moments. See Thorne [ ^6[ and the appendices A and 
B of for many technical formulas and discussion. 

The second step consists of 

(iii) expanding the latter solution in the far wave zone (distance R — > oo) 
in order to find the observable moments of the radiative field which are actually 
measured in the far away detector. 

This step involves constructing a suitable (Bondi-type) coordinate system 
valid in D^, and correcting in particular for the logarithmic deviation of the true 
light cones with respect to the flat light cones. The construction to all orders 
in the post-Minkowskian expansion of such a coordinate system can be found 
in ji^. This yields the operational link between the observable moments and 
the moments parametrizing the multipole decomposition (ii). Computations of 
the observable moments with increasing precision can be found in the sections 
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II.C of H, IV.C of g, V of |6), and in 

Finally it remains to bridge the field in Df. to the inner field of the source 
which is solution in Di of the non-vacuum Einstein field equations. To do so 
we 

(iv) compute the post-Newtonian expansion of the external post-Minkowskian 
field and match it in the overlapping region Di n D^, to an inner solution ob- 
tained by direct post-Newtonian iteration of the field equations in Di . 

This is an application of the method of matched asymptotic expansions [S^ . 
See for a discussion on how useful is this method for bridging the near 
zone and the exterior zone. After matching, the field is determined everywhere 
inside and outside the source as a functional of the source's parameters. This 
allows the determination both of the exterior multipole moments as integrals 
over the source, and of the radiation reaction forces acting inside the source. 
Hatchings have been performed in the sections III of ||5^ , IV o f ]5l| | and III of 
[m] ] to find the multipole moments, and in the sections VI of |49| and in ||5^ 
to find the reaction forces. 



3. PROPERTIES OF RADIATIVE GRAVITATIONAL FIELDS 
3.1. The field in the exterior domain De 

In De we have the vacuum field equations R^i,[g] = 0. It is advantageous to 
use instead of the covariant metric g^i, (with determinant g) the densitized 
metric deviation /i^'' = -y/— gg^" — 77^'' as a basic field variable, where rj'^'^ = 
diag(— 1, 1, 1, 1) is the Minkowski metric. The convention is that the indices on 
h'^'^ (but not on g'^'^ which is the inverse of g^.^) are lowered using the flat metric 
?7/ji> — diag(— 1, 1, 1, 1). With this convention one defines for instance h^,y = 
VfipVi^a-h^'^ and h = rj^^h^'^ . Imposing the harmonic (or De Donder) gauge 
condition for the coordinates a:^ (considered as forming four scalars a;*^, • • • , x^) 
yields = UgX^' = -g^^V^ = -{-gy'^l^d^hf^" , hence 

9^/1^'' = . (7) 

The vacuum field equations in these coordinates then read 

Uh'"' = {h) , (8) 

where □ = □,,= ■q^'^d^di, denotes the Minkowski (flat) d'Alembertian operator, 
and where A^'' (h) can be viewed as a source for the gravitational field in vacuum 
which originates from the gravitational field itself. We have dyK^^ = as a 
consequence of the Bianchi identities, which is implied here by the harmonic 
coordinates condition. 

The source A^^^{h) encompasses all the nonlinearities of the field equations, 
and is at least quadratic in h and its first and second space-time derivatives. For 
instance the quadratic part of A^''(/i), which is such that K^^'^ih) = N^^''{h, h) + 
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0{h^), is explicitly given by {t^^"''^ = + 1"^)) 



^2d^php^dph'''^'' + dah'^pid^K + dph""") 



(9) 



-^dxhp^d^hp'^ + ^dphdph + ^dpKxd'^hp^ 



The cubic part of KP'^{h) can be found e.g. in the section LB of [ p4[ . 

The main technique we use in this section is to look for the solution of the 
field equations in the form of the post-Minkowskian expansion 

hP" ^GhP" + G^hP'' + ■■■ + G" + ■■■ , (10) 

where we denote by G an ordering (or book-keeping) parameter which shall 
be identified later with Newton's constant. The mathematical status of the 
post-Minkowskian expansion is known in vacuum (as is the case here). This 
expansion is asymptotic, in the sense that any solution of the perturbation 
equations (written in (|lll)-(|l^) below) comes from the Taylor expansion when 
G ^ of a family of exact solutions ||8^. Therefore we do not loose any gen- 
erality in assuming the expansion (|l^) , provided that this expansion is carried 
out up to infinity. 

Substituting the asymptotic expansion (|o|) into the vacuum Einstein equa- 
tions (0)-(^, and equalizing all the coefficients of the G"'s in both sides of the 
equations, one obtains an infinite (n = 1, • • • , cxi) set of perturbation equations 
to be solved for all the hP^^s. These equations take the form of wave equations 
in the ordinary sense supplemented by the condition of harmonic coordinates, 

□/ir = Ar(/ii,---,/i„-i) ; 9./ir = 0, (11) 

where the important point is that since A^" is at least quadratic in h (see 
@), all the nonlinear sources A^'' depend for any given n only on the previous 
iterations hl^ with 1 < m < n — 1. For instance we have Aj'' — NP'^ {hi, hi) 
where N^'^ is given by (|9|). 

When n — 1, the corresponding source is zero, A^'' = 0, and therefore 
the coefficient of G in the expansion (|l^), which corresponds to the linearized 
approximation, satisfies 

□/i^ = ; 5^^r = ■ (12) 

An elementary retarded solution in of the source-free wave equation dip = 
reads if = F{t — r/c)/r, where F is an arbitrary function of the retarded time 
t — r/c (with r — |x| the distance to the spatial origin of coordinates located 
in the source). The most general retarded solution is readily obtained by dif- 
ferentiating with respect to space the elementary solution an arbitrary number 
of times, say i : (pi^.-.i^ = dt^ ■ ■ ■dii{F(t — r/c)/r}. We shall extensively use 
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below the notation Ol = di-^ ■ ■ ■ di,, to denote a product of t space derivatives, 
where L ~ 11^2 • • • represents a multi-index composed of t indices. With 
this notation ipi, — dL{F{t — r/c)/r}. Thorne has shown that the most 
general solution h^'^ of the equations ( p^ ) can always be written, modulo an 
infinitesimal gauge transformation preserving the harmonic gauge condition, 
in terms of two and only two types of arbitrary functions F^, which can be 
chosen to be symmetric and tracefree (STF) with respect to the I indices they 
carry (i.e. Fi^ is symmetric under the permutation of any pair of its indices, 
and Si^i^F^ = 0). We shaU denote by ML{t) and SL{t) these two types of STF 
tensorial functions of time and we shall call them the "algorithmic" multipole 
moments of mass type (M^, having ^ = 0, 1, • • •) and of current type (S'l, having 
i = 1,2, • • •). (The adjective "algorithmic" refers to the algorithmic construc- 
tion below of the external field.) In the time-independent case (for stationary 
systems), these multipole moments are identical, aside from normalization, to 
the Geroch and Hansen |89) multipole moments (see the proof in [pO[). 

The moments Mi(t) and SL{t) are for the moment arbitrary except that if 
the harmonic gauge condition in ( p^ ) is to be satisfied, the lowest order mo- 
ments M and Mi (mass monopole and dipole moments) , and Si (current dipole 
moment) are necessarily constant. This corresponds to the laws of conserva- 
tion of the total mass, of the position of the center of mass, and of the total 
angular momentum (the linear momentum which is the time derivative of Mi 
is zero because we consider the case of a source which was at rest in the infinite 
past). Furthermore we shall often choose Mi = by translating the origin of 
the coordinates to the center of mass. The mass monopole M is the total ADM 
mass Q of the source. 

We shall write the result for /i^'' in the form 



cxt 
,4 ^ij 



(13) 
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The ^-dependent coefficients in these expressions are chosen so that Ml (t) and 
SL{t) reduce in the Hmit e ^ corresponding to a Newtonian source to the 
usual expressions of the Newtonian moments of the source |46 . For instance, 



M,j{t)^Q,,{t)+0{e^) 



(15) 



where Qij is the Newtonian mass quadrupole (^ . A complete and self-contained 
derivation of the result (|l^)-(|lj) for the linearized metric can be found in Sec- 
tion 2 of 



3.2. Computing the field nonlinearities in 

Suppose that starting from the linearized metric (p^-(p^ one has succeeded in 
solving the field equations ( p| ) for all the nonlinear coefficients h!^ up to some 
given order n— 1 (l<m<?T. — 1). Each of these coefficients h!^ will be in the 
form of a multipole expansion like h^'^ . How can we find the solution of the 
equations for the next-order coefficient /i^" ? The solution naturally involves 
the inverse retarded d'Alembertian operator defined by 

(□.V)(x, t) = /// ^/(x', ^ - x'l) , (16) 

which satisfies n(n~^^f) = f. However in trying to solve the wave equation 
in ( pi] ) directly with this operator one encounters a problem. Namely the 
source term A^'' is like the previous coefficients h!^ in the form of a multipole 
expansion which is valid only in D^, while the retarded integration in ( p^ 
extends over all space, including the inner domain Di. At the spatial origin 
of the coordinates r located in Di, the multipole expansion is singular 
(see the equations (p^). This problem is of course not a physical problem but 
simply a mathematical one, as we are simply interested in finding a solution of 
the equations (|l|) which is legitimate in Dg. 

As a cure of this problem we use a mathematical trick which allows us to 
find the solution hl^'^ of the equations ( pi] ) in the form of a multipole expan- 
sion, while keeping the range of integration in the retarded integral as it is in 
(p^. This trick consists of replacing the real source A}^'' by a fictituous source 
obtained by multiplying A^" by where i? is a complex number. Assuming 
that the real part of B is large enough (and assuming that the multipole expan- 
sion is actually finite), one "removes" in this way the singularity at the origin 
r — > 0. There is no problem with the behaviour when r — > -|-c» if one assumes 
that the source was stationary in the remote past. Then one can apply to this 
fictituous source the operator (|l6|) as it stands. The result is in general infinite 
at i? = 0, but the point is that the retarded integral of the fictituous source 
has been proved in Section 3 of to admit (after analytic continuation) a 
Laurent expansion near _B — > whose finite part, defined to be the coefficient 
of the zeroth power of B in the expansion, is precisely the solution of the wave 
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equation in (|T^) we are looking for. Denoting the finite part at _B = by the 
symbol FPb=o we pose 

pr = FPB=on]^^(r^Ar) . (17) 

Then pl^ satisfies □p^'' = A^'^ and is in the form of a multipole expansion 
like the previous coefficients /i^'^. However the second equation in (|ll|) which 
is the harmonic gauge condition is still to be satisfied. We shall refer here 
to for the computation from the object pi^^ (or rather from its non-zero 
divergence wi^ = di^pl^'^) of a second object (7^" satisfying at once □g^'' = and 
dvln'^ — ~''^n (^66 Eqs.(4.13) in for the defining expression of q^'^)- With 
these definitions a solution of both the wave equation and the gauge condition 
in (|ll| ) reads 

/^r=pr + 'zr- as) 

This completes the construction of the nth coefficient of the post-Minkowskian 
expansion (px|), and therefore by induction of a whole infinite asymptotic ex- 
pansion solving the vacuum field equations in D^. 

A priori, the above construction seems to be particular, however it has been 
proved in Section 4.2 of |Q to represent in fact the most general solution of 
the field equations in (modulo an arbitrary change of coordinates in D,,). 
The conclusion is that the general radiative field outside an isolated system 
depends on two and only two sets of time- varying multipole moments Af^ (t) 
and SLit) : h'^'^ = h'^'^[ML, Sl]- As a side result but whose proof is in fact 
crucial in the above construction |^^, it is shown that the general structure of 
the post-Newtonian expansion (e 0) of the exterior field in Di n De is of the 
form £'^(lne)P, where k et p are positive integers. This result can be extended 
by a matching argument to the whole near zone Di. Using some energy and 
momenta balance arguments the result could also (probably) be extended to 
the post-Newtonian expansion of the radiative field in The exact level at 
which the first logarithms of e appear is investigated in [|9| (see also |^-[p4|). 

In principle the construction of the field in D^, could be implemented by 
an algebraic computer programme. However more theoretical work should 
be done first on the occurence and properties of high-order tail integrals in 
high-order post-Minkowskian approximations (a preliminary investigation at 
the cubic level can be found in [p7[). 

For practical computations one needs some explicit formulas expressing the 
retarded integral ( p^ when the function / has a definite multipolarity £, say 
/(x, i) — nLf{r,t) where hl denotes the STF part of the product of unit 
vectors til = rii-^ ■■■rii^ {rii = x^/r). As reviewed in and Appendix A 
of the STF tensors = nL{6,4>) can be used equivalently to the usual 
spherical harmonics Y^m{Q,4') for the multipole decomposition. The relevant 
formulas for this are (6.3)-(6.4) in namely 
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where the function g is related to / by 



g{p, s) = I dx ^, 



f{x, s + x/c) 



(20) 



(where a is arbitrary). The formulas (|19[)-(P0[) are correct only when the source 
/(r , t) tends to zero when r — > more rapidly than some power of r (specified in 
|47|). This will always be the case of the fictituous source s/ = r^/ which was 
used in the previous construction of the exterior field (taking Re{B) initially 
large enough and analytically continue the result). Notice the nice correspon- 
dence in (|l9|) between the STF tensor in the source of the wave equation 
and the STF derivative operator 9l (STF part of the product of derivatives 
di^ ■ • ■ di^ ) in the corresponding solution. 

Here are some formulas derived from (|l|)-(|2^) in Section 4 of || and which 
are particularly useful in the study of the quadratic nonlinearities, in which case 
one needs to integrate some source terms of the form f(r, t) = r^^F{t — r/c) 
with k> 2. When fc = 2 we find 



TiLr 



F[t-- 



2(1 



+ CXD 



dXF{t~ X/c) 



(A - rY ln(A - r) - (A + r)^ ln(A + r) 



(21) 



which involves a non-local integral extending over the whole past of the source 
(see below for the associated physical effects). On the contrary the case 3 < 
k < £ + 2 yields a local result : 



'Fit-- 
c 



->k-3 



-3 



{k - 3)\{£ + 2 - k)l 

{e + k~2)\ 

{e + i)\ F'-'^-^-'^t-r/c) 



^ 2Hi(e- 

1=0 ^ 



(22) 



In these two cases the result is in fact finite at i? = and therefore we have 
dropped the Finite Part symbol. When k > £ + 3 the result is no longer finite 
(there is a simple pole at i? = 0) and the expression of the finite part is more 
complicated. See ||5^ for other formulas enabling one to investigate the cubic 
nonlinearities of the field equations. 



3.3. The field in the far wave zone _D„, 

The study of "asymptotics" (namely what is the structure of space-time at large 
distances from the source ?) is an active research field in general relativity since 
the work of Bondi et al [D and Sachs showing that it is possible to find some 
solutions of the field equations with bounded sources in the form of a power 
series in the inverse of the distance R to the source (with a null coordinate u 
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remaining constant). The reformulation of this approach in geometrical terms 
by Penrose |9^ led to the concept of an asymptotically simple spacetime, 
which is by definition a spacetime sharing with Minkowski's spacetime some 
(mathematically precise) global and local asymptotic properties, and which we 
would like to be associated with any radiating isolated system (see [Q-^^ for 
reviews, and the lecture of J. Bicak in this volume). 

The harmonic coordinate system used for the previous construction of the 
field in Df, was very convenient (because the equations for all components of 
the field take the form of wave equations), but it has one drawback : it is not 
of the "Bondi-type" because the expansion of the field when r +00 with 
t — rjc = const involves besides the normal powers of r^^ some logarithms 
of r. This was noticed long ago by Foc k js^. Madore |100[ proved that the 
harmonic-coordinates linearized metric (|l3|)-(p^) cannot be the first approxi- 
mation of a Bondi-type expansion at infinity. The origin of the problem is in 
the deviation of the flat light cones in harmonic coordinates with respect to the 
true null cones. The logarithms of the distance present in 1?^, imply that the 
nth post-Minkowskian approximation there becomes larger than the (n — l)th 
one. However this is not a real problem because the logarithms of r are not 
genuine - they can be removed by a change of coordinates in (in contrast 
with the logarithms of c in the post-Newtonian expansion which are probably 
present in all coordinate systems covering Di). 

The construction of Bondi-type (sometimes also called radiative) coordinates 
from the harmonic coordinates has been achieved in the present frame- 
work to all orders in the post-Minkowskian expansion [Q . A "radiative" post- 
Minkowskian metric was constructed first in the manner of the previous section, 
but correcting at each post-Minkowskian order the flat cones associated with 
the coordinate system in order that they agree (at least asymptotically) with 
the real light cones u — const. Then it was shown (section 4 of jisj) that 
the radiative metric so constructed indeed differs from the general harmonic- 
coordinates metric by a mere coordinate transformation. As a corollary it was 
proved that the Penrose definition of asymptotic simplicity (in a form given 



by Geroch and Horowitz |101]) is formally satisfied for any isolated radiating 
source. Note that this proof holds only in the case where the source is sta- 
tionary in the past, i.e. when there exists a finite initial instant — T at which 
the source starts vibrating and emitting the radiation. If this assumption is 



relaxed, the definition is probably no longer valid 1 102 ]-[ 1041 (at least in its 
original form). 

The main interest of the construction of Bondi-type coordinates in |4^ is that 
it is explicit and can be performed to all orders in G. To leading order in G one 
finds that the retarded time T — R/c'm Bondi-type coordinates (T — R/c and 
the null coordinate u coincide in the limit R 00) differs from the harmonic 
coordinate time t — r/c by 



T-^..-^-^ln 



(^)+0(G^). (23) 
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This is simply the known logarithmic deviation of light cones in harmonic 
coordinates. The constant b is an arbitrary gauge-dependent time scale. The 
correction of order will be given below. 

The Bondi-type coordinates are well-adapted to the study of the field in 
D^, and most importantly permit a clear definition of two sets of multipole 
moments U l (t) and Vl (t) parametrizing the radiative field at infinity, and 
which are directly "measured" in a detector located in D^. We shall refer 
to these moments as the "observable" or "radiative" moments, respectively of 
mass-type (Ul) and of current-type (Vl)- These moments are distinct from 
the algorithmic moments Ml and Sl- Their definition goes as follows. For all 
practical purposes the distance of the source is so large (say R =100 Mpc for 
inspiralling binaries) that one can consider only the leading-order term 1/i? in 
the Bondi expansion. Like in (1) the wave is entirely characterized by the TT 
projection of the spatial metric. It is a standard result (see e.g. |^]) that this 
projection can be decomposed uniquely into multipole moments according to 



1 r 

£=2 

{£ + iyc^-^-^^-b(kVra)bL-2{T - i?/c)| + O (^^^ . (24) 



4G 

/iT.T(X,T) - —V.,kmm NL-2UumL~2{T - R/c 

^ i=2 



2i 



The TT operator Vijkm is given by (||), and t(^km) = ^{tkm + t-mk)- The multi- 
pole decomposition (|2^) is the defining expression for the observable moments 
Ul and Vl (which are functions of the retarded time T — R/c). The total 
power, or luminosity, contained in the waveform (M) is given by 



hoo 



" _^ ..-^"'^ ;(^-i)a!(2^ + i)"'^" ^' 



r2l+l W/' - 1 \PP\(OP + 1)!! ^ ^ 

4^(^ + 2) 



[i - + l)\{2l + l)\\c 



||„2 "L 



(25) 



The ^-dependent coefficients in (24)-(^5|) are ajusted so that the observable 
moments Ul and Vl agree at the linearized order with the £th time derivatives 
of the algorithmic moments Ml and Sl (beware of the fact that Ul and Vl do 
not have the usual dimensions of multipole moments). Comparing (E4| with 
(p^-(p^ (and taking into account a sign difference due to the use in ^13|)-(p^) 
of the metric perturbation \J~qg^^ — 5^^ instead of gij — Sij in (|4|)) we have 

m 

ULiT) = d^ML/dT' + 0{G) , 

Vl{T) = d'SL/dT' + 0{G) , (26) 

where the terms 0{G) symbolize the nonlinear corrections (which can be com- 
puted in principle to all orders using [H). On the other hand we have seen 
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in (15) that the mass quadrupole moment il/y agrees in the Newtonian hmit 
with the quadrupole Qij given by the integral over the source (^. Thus, from 
(ll) and (|l|), we see that in the Newtonian limit the connection between the 
observable moment at infinity Uij and the actual source moment Qij is realized 
(nonlinear corrections in ( p6| ) vanish in the Newtonian limit). Of course this 
yields simply the quadrupole formalism reviewed in Section 2. 

What we shall do in the next subsection is to generalize the relationships 
linking the observable multipole moments to the algorithmic moments in order 
to include explicitly the nonlinear corrections of order G (at least) appearing 
in (p6|). The computation of the algorithmic moments themselves as explicit 
integrals over the source generalizing (15) to post-Newtonian order will be 
dealt with in the section 4. It is clear that knowing both the relations between 
the observable moments Ul,Vl and the moments Ml,Sl, and between the 
moments Ml,Sl themselves and the source, one solves the problem at hand. 



3.4. Nonlinear effects in the radiation field 

It follows from their construction |Q that the observable moments Ul and Vl 
are necessarily given by some nonlinear functionals of the moments Ml and 
Sl (reducing to ( |2^ ) at the linearized order). A very useful information can be 
inferred on dimensional grounds. Indeed the general form of the relations ( p6|) 
can be very easily proved to be 

/~m — l 

n>2 

n>2 ^ 

(endowing the current moment with its natural Levi-Civita symbol), where the 
functions XnL and YnL represent some nonlinear functionals of n moments Ml 
and Sl whose general structure reads 

XnLiT), YnL{T) = 

dvi--- dv„XiT,vi, ■ ■ ■,v„)Mi^^'\vi) ■ ■■Sil"\vr.) . (28) 

J —oo J — oo 

The (a priori) quite complicated kernel X has an index structure made of 
Kronecker and Levi-Civita symbols, and depends only on variables having the 
dimension of time. Accordingly the powers of in each nth nonlinear terms 
of (^ are such that the dimensions of the moments Ul and Vl are respected. 
They turn out (see e.g. the section IV. C of iQ) to be of the type 3{n — 1) -I- 2fc 
where fc is a positive integer representing a number of contractions between the 
indices on the moments composing the term in question. This information is 
useful when determining the type of nonlinear terms which arise in a formula 
like @) below. 
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As expressed in (^8|) the functionals XnL and YnL are given by some highly 
non-local in time (or "hereditary") functionals of the moments Ml and Sl- 
This illustrates the fact that the gravitational field in higher approximations 



depends on all the integrated past (or "history") of the source [105|-|109 
The propagation of the radiation proceeds not only along the characteristic 
surfaces of spacetime or light cones, but also inside them. One can view this 
phenomenon as the propagation at the speed of light of a front wave spreading 
by continuous scattering a secondary wave propagating in average at a speed 
less than c. This wave is referred generally to as the tail of the wave. Perhaps 
a better name would be the wake of the wave. An interesting particular case 
is that of the scattering onto the static spacetime curvature generated by the 
total mass-energy M of the source itself. Often the name of tail of the wave is 
reserved for this specific effect. In that case the tail results from the nonlinear 
interaction between the varying multipole moments describing the radiation 
field, primarily the quadrupole moment, and the static mass monopole moment 
M ' |lO^' H • A signature of the presence of wave tails 



in the gravitational signals received from inspiralling compact binaries should 



be detectable by VIRGO and LIGO [112, 113| 



In order to determine more explicitly the relations ( pTD one must implement 
to nonlinear order the coordinate transformation from harmonic to radiative 
coordinates. This has been done to quadratic order in Section II of To 
see how this works let us first consider the leading-order 1 /r part when r —^ oo 
(with t — r/c fixed) of the linearized metric h'^'^ given by (|l3|)-(|l^). One can 
write 

^r = ;{-^« + -^''(t-r/c,n)|-H0(l) , (29) 

where we have isolated the static term depending on M from the non-static 
contribution z'^" given by 

z-(t,n) = -4E^<W. 

£>2 



i>2 



z-(t,n) = -4^^Mg_,(t) 



l>2 



l>2 ^ >' 



Following the construction of the exterior field in Section 3.1 we insert (29) into 



the quadratically nonlinear source (O) and arrive at the equation to be satisfied 
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by ht^'^ at the leading-order in 1/r, 



□ /if 



1 (AM 



,(2)^1/ 



n - r/c, n) + O 



(31) 



where fc'^ = (Ijii) denotes the outgoing null direction and where 11 is propor- 
tional to the power dEi /dtdQ (per unit of steradian) contained in the linearized 
wave hi : 



1 

4 ^ " Gc^dtdfl 



(32) 



Constructing the solution of ( pl| ) using the method of Section 3.1 necessitates 
knowing the retarded integral of a source with multipolarity £, and behaving 



in r like 1/r^. The relevant formula is (21). Working out this formula to 
leading-order in 1/r (or \nr/r in fact) we get p2| 



-2 



JO 



+ 00 



lyF {t 



In 



(l)+E 



(33) 



showing explicitly how some logarithms of r are generated at the quadratic 
nonlinear level. All these logarithms can be get rid off by implementing the 
coordinate transform of [Esl . With this approximation one finds 



where the vectors and A'^ are given by 



(34) 



\cbJ ' 

t — rj c 



dv n(w, n) 



(35) 



To order G the coordinate transformation (|34|)-(|35|) implies (|23|). Both the 
vectors f ^ and A^ correct for the deviation of the curved light cones from the flat 
cones in harmonic coordinates. After coordinate transformation the metric is of 
the Bondi-type and its 1/R term at infinity is readily obtained. By comparison 
of the result with the expression of the waveform (gj) one obtains the looked- 
for observable moments in terms of the moments , Sl- To quadratic order 
and keeping only the contributions of non-local (or hereditary) integrals one 
finds isl 



Ul{T) = m|'^(t)+ ^3 



2GM 



I ..ln(|)Mr.(T-„ 

/ dvIlLiv) + G1l{T) + 0{G^ 

J —OQ 



2{e +!){£ + 2) J_, 
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+GJl{T) + 0{G^), (36) 

where Hl denotes the STF coefficient in factor of ul in the STF spherical 
harmonics decomposition of 11 (namely 11 = Sn^n^), and where we denote 
by Xl and J7l some instantaneous (by opposition to hereditary) contributions 
which are not controlled at this stage (but see below). 

The expressions (|3^) display two different types of non-local integrals. The 
first type, arising in both Ul and Vl and having M in factor, represents the 
contribution of the backscattering of the waves emitted in the past onto the 
static background curvature associated with the mass of the source (tail effect). 
The constant b in the logarithmic kernel of the tail integrals is the same as in 
(23) and (^. The second type of non-local integral, arising only in Ul (at 
this order), is due to the re-radiation of gravitational waves by the distribution 
of stress-energy of linear waves. We shall refer to this non-local integral as the 
non-linear memory integral. The qualitatively different nature of the tail and 
non-linear memory effects is more clearly understood when taking the limit 
T +00 corresponding to very late times after the source has ceased to emit 
radiation (in the sense that the {£ + l)th time derivatives of the moments 
and Sl tend to zero). In this limit the tail integrals (as well as the instantaneous 
terms Tl and JTl) tend to zero while the memory integral tends to a finite limit 
(hence its name) : 

^ £+1 /?) r~\-oo 

Ul{+^) ^ Mf (+^) + J ^ dvnUv) + 0{G') . (37) 



The nonlinear memory term was first noticed by Payne |114| in a particular 
context. It appeared in the form ( ^6| ) in the present author's habilitation thesis 
(and later in ||5^). Then the memory effect was obtained by Christodoulou 
using asymptotic techniques, and Thorne |117| pointed out that the effect 
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can be recovered simply by adding to his formulas for the linear memory [llS] 
the contribution of gravitons. Wiseman and Will IllE] evaluated the term in 
the quadrupole approximation and found a result equivalent to the ones of 
115: |5^. The nonlinear memory term is however not very important in the 
case of inspiralling compact binaries as compared with the terms associated 
with tails (and even tails of tails). 

It can be readily verified that in a post-Newtonian expansion e — > the tail 
integral present in the observable moment Ul dominates the memory integral 
(see section III. A in 52 ). Indeed the tail arises dominantly at order 1.5PN 



(this is clear from the factor c~ in front of this term), while the nonlinear 
memory term is smaller, being at least of order 2.5PN. For instance, one finds 
that in the post-Newtonian expansion e ^ 0, the mass quadrupole observable 
moment Uij (which is the most important moment to obtain) reads, up to 
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1.5PN order, as 



2GM 



y_ 

.2b 



11 

12 



M[f{T-y) + 0{e'^). (38) 



The constant 11/12 in the integral is actuaUy in factor of an instantaneous 



term (depending only on T through the third time derivative (T)). 
computed in the appendix B of fs^ . 

Pushing the post-Newtonian expansion beyond the 2.5PN order controlled 
in ( ^6|) requires investigating some cubic nonlinearities in the radiation field. 
To 3PN order there appears an effect due to the "tails of tails" , which are the 
tails generated by backscattering of the (quadratically nonlinear) tails them- 
selves onto the background curvature associated with M. This effect corre- 
sponds to the trilinear interaction between two mass monopoles M and (say) 
the quadrupole Mj j describing the linear radiation. The dominant tails of tails 
are computed in |57[ |. Let us give now the expression of the observable mass 
quadrupole moment up to the 3PN order (< ij > denotes the STF projection) : 



It is 



M^f\T) + 2 
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+ 00 
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dv Mi% iv)M^f^, (v) + aMi% {T)Mf^, (T) 



+ mi% {T)M]t^, (T) + jM^,% {T)M,^k [T) + XS^M^:'^^ (r)e,>. 
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2h 



Mlf{T-y) 



(39) 



The various terms are easily recognizable on this expression. Besides the dom- 
inant tail term of order w c~'^, there is the lowest order contribution due to 
the nonlinear memory which is the integral appearing at order (it depends 
quadratically on the quadrupole moment |115, 119, p3), and there is the tail 



of tail term which is given by the last integral of order (one assumes that 
b enters ( p7| ) in the combination {r/cb)^). The a, /9, 7, A, and tr are some 
purely numerical constants in factor of instantaneous terms (like the 11/12), 
and which are computed in . Note that the precise numerical value of these 
constants is not very important physically (they partly reflect our use of har- 
monic coordinates), but their computation using the formulas of Section 3.2 
(and other formulas) is in general long and tedious. 

The observable multipole moments other than the quadrupole admit some 
expressions similar to (|3^), but thanks to the fact that the multipolarity of a 
contribution scales with 1/c in the waveform (24) and with 1/c^ in the energy 
loss (25 ), the needed accuracy for these moments is less than for the quadrupole. 
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We quote here the expressions of the mass octupole Uijk and current quadrupole 
Uij with an accuracy sufficient to include the dominant tails only : 





c Jo 



Mi;>{T~y) + 0{e% 
s\f{T-y) + 0{e% 

(40) 



The computation of the constants 97/60 and 7/6 can be found in Appendix 
C of Q. We assume in that our coordinate system is mass-centered, 
M, = 0. 



4. GENERATION OF GRAVITATIONAL WAVES 
4.1. Post-Newtonian iteration of the field in Di 

We now come to grips with the last step of the present approach, which is to 
find the expressions of the multipole moments Ml and 5*^ as explicit integrals 
over the distribution of stress-energy of the matter fields (and the gravitational 
field) in the source. This problem is evidently crucial in this approach since 
without a satisfactory solution the formulas (^)-(^) giving the observables 
measured far from the source would remain devoided of meaning, and the pre- 
cise interpretation of the gravitational signals received by the future detectors 
would be impossible. 

To be consistent with the formula (^9|) for instance, we should obtain the 
mass quadrupole moment Mij with a precision corresponding to the 3PN order, 
thereby generalizing the Newtonian relation ([l5| ) to take into account all the 
relativistic corrections up to the very high order e^. As we shall see the 3PN 
precision in this moment is not available presently, but My is known with 
the precision which precedes immediately, namely the 2.5PN or precision 
|54[ ^. 



The non-vacuum Einstein field equations in harmonic coordinates read 

d^h'^" = 0, (41) 

where T^^'^ = \g\T^'^ represents the stress-energy distribution of the matter 
fields (with T'^^ the usual stress-energy tensor), and where the other notations 
are defined in Section 3.1. In the near-zone Di the field equations can be solved 
by means of a direct post-Newtonian iteration when e ^ 0. We shall perform 
it without making any assumption concerning the stress-energy distribution 
T^'' (it could be for instance that of a perfect fluid, or it could describe point 



HOUCHES 



23 



a = — 

c 



particles without interactions), except that it should have a compact support 
and correspond to a source whose parameter e is of the order 0.3 (say) at most. 

It is useful to introduce as basic variables describing the source some densities 
of mass fj, current Oi and stresses defined in terms of the contravariant 
components of the stress-energy tensor T^^ (with T" = SfJ^-T'^ ) : 

c 

= T'' ■ (42) 

In particular a is related to the Tolman mass density valid for stationary sys- 
tems Q|. Here the powers of are such that these densities have a finite 
non-zero limit when e — > (T^"^ has the dimension of an energy density) . The 
usefulness of these definitions lies in the fact that they simplify appreciably the 
first post-Newtonian approximation (this was pointed out in Section II of [^), 
and therefore, as we shall see, the subsequent approximations built on it. See 
the lecture of T. Damour in this volume for the derivation of the IPN approx- 
imation using the source densities (^2|), and for applications to the relativistic 
N-body problem. 

From the densities (^) we introduce next the following retarded potentials 

V = -47rGn^V, 
V, = -47rGn^V,;, 



^ Uvd.v - ^s,,dkVdkV 



(43) 



where is the flat retarded d'Alembertian operator (p^). The scalar and 
vectorial potentials V and Vi reduce in the limit £ — + to the Newtonian 
potential U and the usual "gravitomagnetic" potential Ui. The tensorial po- 
tential Wij is more complicated but is generated in the limit e — > by the 
sum of the matter stress density in the source and the stresses associated with 
the Newtonian gravitational field itself (the gravitational stresses are of the 
same order as the matter stresses when e ^ 0). The sources of the potentials 
V and Vi are of compact support, while the source of Wij is non compact. 
Thanks to the equations of continuity and of motion for the source densities 



(42) these potentials satisfy some simple differential identities to Newtonian 



order : dtV + d^V^ = 0{e^) and dtVi + djW^j = O(e^). Beware of the fact 
that these potentials are different from the external potentials introduced in 



(li3|)-(y). 

To Newtonian order one has = -AV/c^ + 0(£^), h"" = 0{e^) and h'^ = 
Oje'^), that we can substitute into the right-hand-side of the wave equation in 



(41 ), where with this approximation the nonlinear term K^'^(h) can be replaced 
by its quadratic nonlinear piece given by (^. Applying to the result 
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(without any finite part of course because the source is regular in Di), and 
neglecting consistently all the small terms, one obtains the first post-Newtonian 
approximation in the form 

C = -^V+^iWu'2V') + Oie<^) , 

= -^W.,+0{s'). (44) 

Actually this approximation is slightly more accurate than the IPN approx- 
imation because it takes into account the terms of order not only in the 
00 component of the metric, but also in its ij components. We indicate by 
the subscript "in" that the metric is the inner metric valid in Di, which dif- 
fers by a coordinate transformation from the exterior metric of Section 3 (see 
below). Substituting the expressions (^J) back into the right-hand-side of the 
field equation in ([4l|), where now the nonlinear term K^^{h) should include 
besides the quadratically nonlinear piece (^) a cubically nonlinear piece (given 
in Section LB of [Q), neglecting consistently the higher order terms and then 
inverting the result by means of D]^^, one gets 



0(e^£^£«), (45) 



where the remainder means that the neglected terms are O(e^), O(e^) and 
0{£^) in the components /iP°, /iPj and /i-;^, respectively. The quantities T^'^ and 
are given by the corresponding quantities in the right-hand-side of 
(pl|), but where one retains only the terms in the post-Newtonian expansion 
up to the accuracy of the remainder in (|4|). These are given as explicit com- 
binations of derivatives of V , Vi and Wij. For instance the effective nonlinear 
source A*^ (F, W) reads 

A°°(y,W) = -^dkVdkV+^i^-Vd^tV-2VkdtdkV 

-Wkmdl^V + ^(dtV)^ + IdkVraidkVm + 
o Z 

+dkVdtVk + 2dkVdkWram " IvdkVdkV 



A°'(T/,w^) = ^i^dkV{d,Vk-dkVi) + -^dtVd,v^ , 

IC' (1/, w) ^ ^ \d,vdjV - ^s,,dkVdkV^ + |2%yatFj) - d.VkdjVk 
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+ 7i5,jdMdkV^ - dM \ ■ (46) 



See the section II. A of [p4| for the derivation of these formulas. 



4.2. The multipole moments as integrals over the source 

The expressions of the multipole moments Ml and Sl will follow from the 
comparison between the inner metric (^5[)-(p6|) valid in Di and the re-expansion 
when e ^ of the exterior metric valid in D,,. The comparison (and in fact 
matching) is done in the overlapping region D^. n Di which always exists for 
slowly-moving sources. This "matching" region is where one can simultaneously 
expand the field when e — > and decompose it into multipole moments. The 
matching requirement is simply that the two inner and outer post-Newtonian 
asymptotic expansions should be (term by term) identical in Di n D^, after 
the performing of a suitable coordinate transformation. Let this coordinate 
transformation be 

^xt^-^ + ^^l^in) , (47) 

where x^^^ and denote the coordinates which were used in D^ and Di 
respectively (with a slight change of notation with respect to Section 3). 

Recall that the exterior metric is constructed starting with the linearized met- 
ric hf^^ whose components are defined by some external potentials V'^^^ , Vi^^ 
and V^™* given explicitly as some infinite multipole expansions parametrized by 



Ml and Sl (see (|13|)-(14)). The problem is to relate via the matching require- 
ment (^^ these exterior potentials to the inner potentials V , Vi and Wij used 
in Di. In fact we do so first with some intermediate accuracy. The method is to 
look for a numerical equality in Di n Df. between exterior and inner potentials, 
and then transform this equality into a matching equation (i.e. an equation re- 
lating two mathematical expressions of the same nature) by replacing the inner 
potentials by their multipole expansions in Df,- The matching equation is valid 
formally "everywhere" , and it can be used to get the functional relationships 
linking the multipole moments to the source. 

For instance it is found that the matching equations relating the potentials 
V'^^^ and Vj™* to the multipole expansions M{V) and AA{Vi) of the compact- 
support inner potentials V and Vi (where M refers to the multipole expansion) 
read as 

F^'^t = M{V) + cdtv"" + 0{e^) , 

y-t = M{Vi)~'^d,^'' + 0{e') , (48) 
where ip^ denotes the zero component of the gauge transformation vector in 



( |47[) . The equation for V is valid to first post-Newtonian order, but the 
equation for Vi^'^ is valid only to Newtonian order (it can be proved that is 
of order e^). The matching equations ( ^8|) were obtained in Section II of jscfl . 
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The multipole expansion of a retarded potential with compact- support source 



is known from the work of CampbeU, Macek and Morgan |120|. It has been 



re-calculated in the appendix B of ||50| using STF spherical harmonics. Using 



this formula Damour and Iyer |Q improved much the work |120 dealing with 
the case of linearized gravity. The result (with obvious notation for the index 
H) is 



It/ r 



r " \ c 



(49) 

where the multipole moments which parametrize the expansion are given by 
V!:{t) = J d^^XL J dzSi{z)af,{x,t + z\K\/c) . (50) 

We denote by the STF part of the product of £ spatial positions xl — 
■ ■ ■ Xif. The function 6i{z) takes into account the physical delays due to the 
propagation of the waves inside the compact-support source. It is given by 

'^K^) = ^^TT§^(1-^')' ; f^dz5,{z)^l (51) 

{5^ can be simply related via the Rodrigues formula to the Legendre polynomial 
of order tj. With the two matching equations (|8|) and the explicit formulas 
(pO|)-(^, it is a simple question of algebraic manipulation (using STF tech- 
niques), and analytic post-Newtonian expansion (using the formula (B.14) of 
[50|), to deduce the multipole moments Ml and Sl entering the left-hand-sides 
of (|4^). In particular the mass moment Ml is obtained with post-Newtonian 
accuracy as 

M.(,) ./.3.{i,,+^J^a?.-^|^±|fi±^a.,. j+o(.') ^ (52) 

This expression generalizes to IPN order the result (^5|) of the Newtonian 
formalism. The striking feature about this expression is that the integrand 
has a compact support, althought it includes the contribution of the gravita- 
tional field to IPN order, which is expected to extend over the whole three- 
dimensional space. However, the gravitational contribution to IPN order turns 
out to be entirely contained into the mass density a. As the expression ( p2| ) has 
been derived rigorously within the post-Newtonian framework, it can also be 
said (when combined with the previous results in D^,) to bring a satisfactory 
proof of the lowest-order quadrupole formula itself. 

To the next order (2PN) things are more delicate precisely because of the 
explicit contributions of the gravitational field which come into play, making 
the expressions of the mass multipole moments a priori non compact. This 
problem has been solved in where it was found necessary to obtain first 
the expression of the multipole expansion M {Wij ) of the inner potential Wij 
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defined in (p3|), and whose source we recall is non-compact. Thus the formulas 
(50)-(|T]) cannot be employed in this case, but instead we have obtained (in 
Section III.C of Q) 



e>o 







>5 ( 









( -dM{V)d,MiV) + -S,,dkM{V)dkMiV) 



(53) 



where the multipole moments are given by 
W,5(i) =FPb=o J d'x|x| 



X / dz6£{z) 



XL 

1 

4^ 



^^V^JV - -S,,dkVdkV 



(x,t + z|x|/c) . 

(54) 



The first term in (|5^) together with the moments ( p4D look like exactly what 
we would obtain by using blindly the formulas (50)-(^l[), even though the po- 
tential is non compact. However there is an important difference which is the 
presence in the moments ( p^ of the finite part procedure which was used in 
the definition of the exterior metric (see (p^)). The role of this finite part is 
to deal with the infinite support of the integral ( p4[ ) , which would be divergent 
otherwise (because of the presence in the integrand oi behaving like |x|^ at 
spatial infinity). Notice that the finite part is conveyed all the way through 
the analysis, and is seen a posteriori to be in fact necessary to ensure the con- 
vergence of the integrals. Thus no ad hoc prescription has been used to obtain 
the moments ( |5^ ) in a well-defined form mathematically, even in the case of a 
non-compact supported source (see the proof in ||5^). As for the second term 
in (I53I), it ensures as is also seen a posteriori that Ai{Wij) satisfies the cor- 
rect wave equation (deduced from ([43|)) outside the source. This second term 
involves the multipole expansion M(V) of V in conformity with the retarded 
integral operator FP B=QO]j^r^ which is defined only when acting on multipole 
expansions. 

With this sub-problem of obtaining A4{Wij) solved we can write down sim- 



ilarly to (48) a matching equation relating the external potential V^'j^^ to the 



multipole expansion M{Wij). Together with (|48| ) this permits us to relate the 
nonlinear source in the exterior field, which to this post-Newtonian order is 
shown to be exactly A'^'^ {V^*' , W^^*') , where A'^" is given by (^) and where 
W°j^^ is defined from V°j'^^, to the multipole expansion of the nonlinear source, 

that is A {M{V),M{W)), modulo the terms associated with the coordinate 
transformation. This relation being established one can write a matching equa- 
tion valid to higher post-Newtonian order. The nonlinear part of the coordinate 
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transformation is found to cancel out, so that it remains simply a linear gauge 
transformation associated with the vector ip'^. Finally a reasoning similar to 
the one followed to obtain (p3[)-(p^ leads to the linear piece ft,^"^ of the external 
field (which is a functional of the algorithmic moments Ml and Sl) in the form 
of the explicit multipole expansion 

+ d^'"' + O(e^) , (55) 

where d(p^'^ ~ d^(p'^ + d'^(p^ — ■q^^'^dxip'^ is the linear gauge transformation, 
which depends on some (reducible) multipole moments given by 

T^''{t)=FPB=o J d'x\xfxL J dz(5,(z)T^''(x,t + z|x|/c) . (56) 

The integrand of these multipole moments is nothing but the total stress-energy 
momentum of the matter fields and the gravitational field (with 2PN accuracy) , 
namely 

T^''^T'"' + ^T''iV,W), (57) 

which is given, thanks (in particular) to ( [l6|), a s a fully explicit functional of 
the mner potentials V, Vi and Wij. Like in (|54|) there is a Finite Part symbol 
in front of the integral which is seen (a posteriori) to make the non-compact- 
support multipole moments to be mathematically well-defined. 

The result (p5|)-(p7|) (proved in Section III.D of Q) is satisfying because 
one recovers finally the expression of the total pseudotensor of the matter and 
gravitational fields in the source, which is exactly the one which appears in 
the right-hand-side of the Einstein field equation (^l]), except that here it is 
expanded when e — > up to 2PN order. Notice that although the range of 
integration extends up to infinity, it is the post-Newtonian expansion valid in 
Di of the pseudotensor which is to be used in (|5^). This is thanks to the finite 
part procedure and the properties of analytic continuation (see js^ for the 
reason why this works). 

The form of the result (p5|)-(p7|) is therefore identical to what we would 
obtain by writing from the field equation ( |4l| ) the solution directly as the re- 
tarded integral of the total pseudotensor (this operation is perfectly licit), 
but then by writing uncorrectly the multipole expansion of that solution apply- 
ing the formulas (^)-(|5l]). This second operation is illicit because the formulas 
apply only to compact-supported sources. As a result the multipole moments 
would be given by divergent integrals. This was the approach initially followed 
by Epstein and Wagoner |38| to IPN order, and then generalized formally by 
Thorne to all post-Newtonian orders. The error is not too severe to IPN 
order because one actually recovers the correct result by discarding the infinite 
surface terms in the Epstein- Wagoner moments. This has been proved using 
the compact-support formula ( p2|) (see the appendix A of |Q), and this has 



4G 



Ghr[ML,SL] = -—J2 
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also been checked in applications to binary systems [ p9[ . However in higher- 
orders the neglect of infinite surface terms in the Epstein- Wagoner moments 
becomes physically uncorrect, because it is then equivalent to the neglect of 
physical effects in the waveform, namely the tail and other nonlinear effects 
computed in Section 3.4. This is clear from the fact that by (p5[)-(p7|) one sees 
that the multipole moments which are given by (the finite part of) integrals 
over the total stress-energy pseudotensor are simply the moments and S^, 
which differ from the observable moments Ul and Vl in the wave field (|2j) by 
all the nonlinear effects in D^. 

Recently Will and Wiseman |Q have succeeded in curing the defects of the 
Epstein- Wagoner approach. As a result they obtain a manifestly convergent 
and finite procedure for calculating the gravitational radiation from isolated 
systems to post-Newtonian order. 

It still remains to compute from (|55|)-(|5^) the (irreducible) moments and 
Sl , which means to decompose in the right-hand-side of ( p5| ) the moments JFj^" 
into irreducible STF pieces. Apart from the Finite Part symbol in front of the 
moments and the post-Newtonian remainder in (|55|), this technical problem is 
the same as in the case of linearized gravity in which case it has been solved 
by Damour and Iyer The mass multipole moment Ml to 2PN order is 

therefore deduced, and, after further transformations done in Section IV. B of 
]54| ], we end up with 



ML{t) = 



o- + -^{ciiU - aPii) 



2c2(2<' + 3) 
1 



q2 

* " c2(^+l)(2^ + 3)~ 



4C7\ 



„ 2 , dkU[d^Uk - dkU^] + -dtUd^U 
ttGc^ \ 4 



8c4(2£ + 3)(2£ + 5) 

c^{e+l){i + 2)i2£ + 5) '''' ' 
1 



c4(^-|-l)(2^-K3)(2^ + 5) 



7i3 + -^diUdjU 
■' AttG 



XL 



Pjdf^U - 2U^^t^^U + 2^^UJ^,U^ 



+ 0(£^) 



(58) 



The potentials [/, Ui and Pij are the Newtonian-like potentials to which tend 
the potentials V ^ Vi and Wij in the limit £ — > 0. As we see there are many ex- 
plicit contributions to 2PN order of the non-compact-supported distribution of 
the gravitational field. The expression ( |5^ ) has been generalized recently to in- 
clude the next 2.5PN correction [EgI. Similarly one deduces also the expression 
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of the current multipole moment Sl, but with IPN accuracy only : 

SUt) = FPB=oeab<^, J d3x|x|^|xi_i>, (^1 + Ic/^ a, + ^^^^^iTif ''" 



1 ^ 



{2e+l)XL-l>ac 



dkUidbUk - dkUb) + ^dtUdbU 



dt 



cTbc + -^dbUdJJ 
47rG 



+ 0{e^) . (59) 



The first correct expression for the IPN current muhipolc moment was in fact 
obtained earUer by Damour and Iyer making extensive use of distribu- 
tional kernels to deal with the quadratic nonlinearities of the gravitational 
field. Their result has the advantage over ( |59| ) to explicitly express the IPN 
current moments Sl in terms of compact-support integrals. Its extension to 
higher orders would require to generalize the quadratic kernels they introduce 
to cubic and higher kernels. [In a sense this generalization has been performed 
in Section III.C of where the result (58) has been transformed in an ex- 
plicitly compact-support integral.] See Section IV. B of js^ for the complete 
equivalence of with the result of ISl]]. 



5. GRAVITATIONAL RADIATION REACTION EFFECTS 

On physical grounds one expects that it should be possible to any post-Newton- 
ian order to interpret the results concerning the outgoing radiation field in 
(and the energy, linear momentum and angular momentum it carries) as due 
to some radiation reaction force acting locally inside the source (in Di). Indeed 



this is known to be correct at the Newtonian level |16 , [M-|30[|. In Section 2.2 



is reviewed the Burke and Thorne |20-|^ expression for the reaction force 



at this level. Other expressions of the force in other coordinate systems were 
obtained by several authors, notably Chandrasekhar and Esposito [n6[ (see 



85 for a discussion). Recall that the (Newtonian) radiation reaction force is 
directly responsible for the observed acceleration of the orbital motion of the 
binary pulsar PSR 1913-^16 |5|. 

Unfortunately the problem of radiation reaction onto the source is not so 
well understood as the problem of the generation of radiation in the regions far 
from the source (sections 3-4). Here we shall report the results of the extension 
to first post-Newtonian order (and even 1.5PN order) of the lowest-order Burke 
and Thorne radiation reaction force (6) . These results confirm to this order the 
physical expectation that all effects of the reaction force are contained in the 
fluxes of energy and momenta carried out by the waves at infinity, but clearly 
some more work should be done on this problem. 

The novel feature when one goes from the Newtonian reaction force (6) to 
the first post-Newtonian one is that the reaction potential is no longer com- 
posed of a single scalar (depending on the mass- type multipole moments), but 
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involves also a vectorial component depending on the variations of the current- 
type multipole moments. The vectorial component of the reaction could be 
important in some astrophysical situations like rotating neutron stars undergo- 
ing gravitational instabilities. It was first noticed in the physically restricting 
case where the dominant quadrupolar radiation from the source is suppressed 

(m). 

From the general case investigated in ||5^, ^ it follows in particular that the 
expressions of the scalar and vectorial reactive potentials which generalize to 
IPN order the Burke-Thorne potential (6) are given by 



l/''°^'=(x,t) 



4 

45' 



70' 



(60) 



One recognizes in the first term of the Burke-Thorne scalar potential 

for the mass quadrupole moment My . The vector potential ^^oac depends 
notably on the current quadrupole moment Sij. It is crucial that Mij in the 
Burke-Thorne term should be given with relative IPN precision in order to be 
consistent with the accuracy of the approximation. The correct expression for 
that moment is above (taken for £ = 2). The other moments in (|6^) take 
their usual Newtonian forms at this approximation. Notice that the expressions 
(60) are valid within the full nonlinear theory. They are shown in |^ to arise 
from a specific time-asymmetric component of the nonlinear gravitational field 
in De- The matching of these expressions to the field of the source in Di is 
done in 

The energy, linear momentum and angular momentum which are extracted 
from the source by the IPN radiation reaction force associated with the scalar 
and vectorial potentials (^) are to be computed locally (in Di) using the 
equations of motion of the source. This computation has been done recently 
[5^ with a result in perfect agreement with the corresponding IPN fiuxes of 
energy and momenta already known from their computations in the far zone 
[|6[|50). Namely, 
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(61) 



where E, Pi and Si denote some local (instantaneous) quantities which are 
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given by some compact-support integrals over the matter in the source, and 
agree to IPN order with the usual notions of energy, linear momentum and 
angular momentum of the source. These formulas show the validity to IPN 
order of the (generally postulated) balance equations between the decreases of 
energy and momenta in the source and the corresponding fluxes at infinity. Of 
particular interest (notably because this is a purely IPN effect) is the decrease 
of the linear momentum dPi/dt, which corresponds to a net recoil of the center 
of mass of the source in reaction to the emission of waves. Numerous authors 
40 , 122, 123, 4^ had computed before this effect as a flux of linear momentum 



at infinity. 

In the case of binary systems, Iyer and Will |124, 125| have obtained, assum- 
ing the validity of the balance equations for energy and angular momentum, 
the radiation reaction force in the equations of motion of the binary to IPN 
relative order (which corresponds to the 3.5PN order beyond the Newtonian 
acceleration because the reaction effects arise at 2.5PN order). Their result is 
valid for a large class of coordinate systems. They also show |125[ that the 
reactive potentials ( |60| ) when specialized to binary systems correspond in their 
formalism to a unique and consistent choice of a coordinate system (namely, 
a "Burke- Thorne-extended" coordinate system). This represents a non-trivial 
check of the validity of the reactive potentials (|60|). 

To the next 1.5PN order in the reactive potentials there appears the same 
phenomenon as in the radiation field in D^u, namely the occurrence of hered- 
itary effects associated with tails. This is required if the balance equation for 
the energy is to remain valid at this order. That this is indeed the case has 
been proved in p9[ | where the dominant hereditary contribution in the post- 
Newtonian expansion of the equations of motion of the source (in Di) was 
shown to arise at the 4PN order beyond the Newtonian acceleration, that is 
precisely at the 1.5PN relative order in the reaction force. Therefore, one can 
say that locality in time is lost at the 4PN order (or order e^) in the equa- 
tions of motion of a self-gravitating source. This hereditary contribution can 
be shown in a particular gauge to modify the scalar reactive potential in 
by the tail integral |49| 



+ CXD 



dy 



In 



f-U- 
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(62) 



This tail contribution in the reaction is recovered in Section III.D of ||5^ using 
a more systematic method, and the constant 11/12 is computed. 

The following picture can be given to this phenomenon. The waves emitted 
by the source at all epochs in the past are scattered onto the spacetime cur- 
vature generated by the mass M. The secondary waves which are produced 
converge back onto the source at our epoch and modify its present dynamics. 
The total reaction force is thus the vectorial sum of the reaction force due to 
the present emitted radiation (this is (|o|)), and of the reaction force due to 
the incoming tail radiation (see (|62|)). It is shown in Section III.D of (52j that 
the tail-induced reactive potential SV'"^'^'^ is consistent with the tail integral in 
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(38) in the sense that it gives rise to a term of order in the energy balance 
equation in (^) which is in perfect agreement with the contribution of the tail 
integral to the flux of energy in 



6. APPLICATION TO INSPIRALLING COMPACT BINARIES 

Inspiralling compact binaries are very relativistic (but not fully relativistic) 
binary systems of compact objects (neutron stars or black holes) in their late 
phase of evolution which precedes immediately the final coalescence, during 
which the two objects spiral very rapidly around each other under the influence 
of the gravitational radiation reaction forces. 

It has been recognized in recent years that a very precise prediction from 
general relativity, taking into account many relativistic (or post-Newtonian) 
corrections, will be necessary in order that all the potential information con- 
tained in the signals from inspiralling binaries can be deciphered by the future 



detectors LIGO and VIRGO |^-|l3§,ll]J |l]J. The demanded precision 
can be calculated using black-hole perturbation techniques in the special case 
where the mass of one object is very small as compared with the other mass 



|134| -|140|. This precision turns out to be at least 3PN beyond the quadrupole 
radiation | 136| , 137, 140 1. Thus, inspiralling compact binaries offer us the ideal 



application of the previous relativistic formalism — and constitute presently 
our main motivation for its development. 

The information contained in the signals of inspiralling compact binaries will 
be extracted from the noisy output of the detector by a technique of matched 
filtering, necessitating the a priori knowledge of the signal. Crucial to a suc- 
cessful data analysis will be the knowledge of the instantaneous phase of the 
binary, namely the angle between the separation of the two bodies and a fixed 
direction in the orbital plane. 

The time evolution of the orbital phase should in principle be determined 
from the gravitational radiation reaction forces acting locally on the orbit. 



However, as these forces are only known to IPN order |124, 125, ^ |5^ and 
1.5PN order p9|], one relies in practice on an energy balance equation similar 
to the one in (pTl) , namely 

where in both sides one uses some very precise expressions of the (binding) 
energy E of the binary and its total luminosity C in gravitational waves (the 
angular momentum balance equation is not necessary because the orbit has 
been circularized by the radiation reaction forces). The fact that one is obliged 
presently to postulate the validity of the energy balance equation (|6^) to higher 
order than 1.5PN (for instance 2.5PN as postulated below) is certainly a weak 
point in the analysis, which shall have to be improved in future work. 

The binding energy E of the binary which appears in the left-hand-side 
of ( |63| ) results from the equations of motion of the binary which have been 
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obtained to 2.5PN order by Damour and Deruelle |3^-|^ in their study of 
the binary pulsar PSR 1913+16. For inspiraUing compact binaries one simply 
needs to specialize these equations to the case where the orbit is circular (apart 
from the gradual inspiral). The 2.5PN-accurate binding energy reads [ p2| 

E = -^Mi^xi^l-^ ^(9 + v)x-^(27 -19iy+Y^ + Oix^)'j . (64) 

(See also the appendix B of for a derivation for this formula from the 
Damour-Deruelle equations of motion.) The total mass of the binary is M = 
mi + m2, and v = mim2 / M'^ represents a particular mass ratio which is such 
that < V < 1/4, with i' — 1/A corresponding to the case of two equal masses 
and 1/ ^ corresponding to the test-mass limit for one body. For convenience 
we denote by a: = {GMlj / c'Y^^ a small dimensionless post-Newtonian param- 
eter of order e^, where to is the instantaneous angular frequency of the orbit 
(uj = Itt/P, where P is the orbital period). 

On the other hand, the 2.5PN~accurate luminosity C of the binary is obtained 
by application of the formalism above. Namely, one uses the decomposition of 
C in terms of observable moments Ul and Vl as given by (p5|). Then one 
relates the observable moments to the moments and S'l by the equations 
(^6|)-(^0|). Actually to 2.5PN order the mass quadrupole moment comprises the 
dominant tail integral and the quadrupole-quadrupole interaction terms which 
include the nonlinear memory integral, but the tail of tail integral appearing 
in ( ^9|) is negligible. Nonlinear effects are rather easy to compute for binary 
systems thanks to the study done in |Q , where the required formula is proved 
to apply to a binary which is actually spiralling (or decaying) . Indeed this is not 
obvious a priori because we are considering hereditary effects which depend on 
the whole integrated past of the binary (see the appendix B in [|o)). In a second 
stage one must use the expressions of the moments Ml and Sl as integrals over 
the source. These are given by (^) and ( ^9| ) (to 2.5PN order the next-order 
term in Ml is a priori also needed, but in fact this term is instantaneous and 
does not contribute in the energy loss to this order, see below). A hard task, 
carried out to 2PN order in and |6^, is to reduce all these integrals for 
compact binary systems. This entails using delta functions to describe the 
compact objects (but this can probably be justified using the results of pj- 
35 ). Notably the main cubically nonlinear term in ( p8| ) (term with Pijdf^) 



deserves a special treatment. This term has been computed by two different 
methods in |Q (see Section III.C and Appendix B there). The final result for 
£ reads 



32c5 2 sf, /1247 35 \ ^ 3/2 
C = — —v x^il~i \ ly ] X + Attx-^'^ 



5G 1 V 336 12 

44711 9271 65 , , 2 

+ 1 V -I v'' X 

' 9072 504 18 ' 
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/8191 535 ,^^5/2 ^^^^3 ^g^^ 



V 672 24 

This expression was already known to Newtonian order by Landau and Lif- 



chitz |8C|, |8l|], and to IPN order by Wagoner and Will |141| using the Epstein- 
Wagoner moments p8[ (in which infinite surface terms are to be discarded). 
It has been re-calculated to IPN order in |^ using the well-defined compact- 
support moment (|5^). The 1.5PN order is the 47r term which is the contribution 
of the dominant tail integral present in (|3^). The coefficient 47r is specific to 
binary systems moving in circular orbits, and was obtained first by Poisson 
]135| using a black-hole perturbation technique valid in the test-mass limit 
V ^ (this happens to give the answer, which does not depend on even 
for arbitrary mass ratios). The computation of the tail integral in ( ^8|) in the 
case of binary systems was done in j6^, Q with result yielding the 47r term. 
The 2PN term computed in |]6l| is entirely due to relativistic corrections in the 
moments of the source, most importantly the mass quadrupole moment ( ^8|) 
(issued from |Q). Will and Wiseman computed independently this term 
using the Epstein- Wagoner moments which give the correct result, as proved 
in |5^, because there are no tails at this order (provided that infinite surface 
terms are discarded). (But see [|4| for improvements of the Epstein- Wagoner 
approach, and for how to compute tails within this approach.) The common 
result of these 2PN computations was summarized in |62| , where it was also 
pointed out that the effects of intrinsic rotation of the two objects are in general 
negligible as compared to the gravitational effects. Finally the 2.5PN term has 
been added in [56| and shown to be due exclusively to the mass octupole and 
current quadrupole tail integrals in ( ^ ) , and to the post-Newtonian corrections 
in the dominant mass quadrupole tail integral in (|3|) or (p9|). The fact that 
only the tail integrals give a contribution in C to 2.5PN order is special to the 
case of a circular orbit [ p6| . In the non-circular case the instantaneous terms 
would also contribute. Incidentally notice that the nonlinear memory integral 
is instantaneous in the energy loss, and therefore gives no contribution in (|65|). 
Thus the nonlinear memory effect has rather poor observational consequences 
for inspiralling binaries. 

The tail of tail integral given by the last term in (^9|) gives a contribution 
to £ at the first order which is neglected in (|65|), namely the order 3PN or x^. 
This contribution to C is computed in ||57| (it involves a logarithmic term in the 
frequency) . However the whole 3PN term in ( |65| ) will include also many contri- 
butions coming from the relativistic corrections in the moments of the source. 
These contributions are not under control presently. The next 3.5PN order is 
expected like the 1.5PN and 2.5PN orders to be due entirely to higher-order 
multipolar tails, and to post-Newtonian corrections in low-order multipolar 
tails. The still higher 4PN order includes many effects and seems presently 
difficult to reach within the post-Newtonian theory. 

There is a case, however, where C is known up to 4PN order. This is the 
test-mass limit case v ^ 0, which has been computed by black-hole pertur- 
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bation techniques. Tagoshi and Nakamura |137| first obtained numerically the 



coefficients in C up to 4PN order. Then Sasaki |138| showed how to solve 



in an iterative way the Regge- Wheeler equation for the perturbations of the 



Schwarzschild black-hole, and finally, based on this, Tagoshi and Sasaki |139| 



obtained analytically the coefficients in C up to 4PN order. (See also [142, 143| 
for the case of a rotating black hole.) The result ( |65| ) of post-Newtonian theory 
is in perfect agreement in the limit v with the truncation to 2.5PN order of 
the result of perturbation theory (tails of tails computed in also perfectly 
agree). As the post- Newtonian and perturbation theories are so different, the 
agreement is very satisfying, although mandatory of course. 

Finally it remains to substitute the 2.5PN expressions of both E and £ into 
the energy balance equation (|63|) (assumed to be true at the 2.5PN level), 
in order to obtain the instantaneous phase of the binary as a function of 
time. The phase is = J ojdt where w is the orbital frequency, from which 
one deduces = — / {lo/ C)dE. The result is simpler if we use instead of the 
local time t flowing in the observer's frame {t = T — R/c with the notation of 
Section 3) the adimcnsional time variable [ |6l| 

e^^(t.-0, (66) 

where tc denotes the instant of coalescence (at which the frequency goes for- 
mally to infinity). In terms of this time variable the orbital phase is obtained 
as 



V8064 96 y 4 
/ 9275495 284875 1855 



V 14450688 258048 2048 
/ 38645 15 



V Trine -HO(e"i/«) } , (67) 
V172032 2048 J ^ ;j > v ; 

where (po is some constant phase determined by initial conditions (for instance 
when the signal enters the frequency bandwidth of the detector). This expres- 
sion is valid only in the post-Newtonian regime where = 0{e^). 

Note that the determination of the phase (|6^) does not constitute by itself 
the complete answer to the problem because the waveform (described by two 
independent polarization states h+ and h^ ) is also to be computed. Thus the 
phase (/)(t), and the orbital frequency uj{t) = <j){t), have to be inserted into 
the waveform of the binary, which is itself to be known with the best possible 
post-Newtonian precision (though the latter precision does not need to be so 
high as in the determination of the phase; see, e.g., |131[ |). It is known Q that 
the waveform introduces in the phase evolution some contributions which add 
to those coming purely from the energy balance equation (|6^) , but these extra 
contributions arise at an order which is higher than 2.5PN. For the computation 



of the waveform to 2PN order see |63 
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